This program calculates probabilities for observing an event at least once, given the true probability inputs. Creates a table that can display multiple probabilities.
The program is written in R.
View Code
probability_of_observing_a_rare_event = function(event_rates, sample_size)
{
rates = sort(unique(event_rates))
probabilities = sapply(rates, function(event_rate) {
paste0(round((1 - ((1 - event_rate)^sample_size)) * 100, 1), "%")
})
result = list(
n = rep(sample_size, length(rates)),
rate = rates,
probability = probabilities
)
return(result)
}